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We prove that the results of a finite set of general quantum measurements on an arbitrary di- 
mensional quantum system can be simulated using a polynomial (in measurements) number of 
hidden-variable states. In the limit of infinitely many measurements, our method gives models with 
the minimal number of hidden-variable states, which scales linearly with the number of measure- 
ments. These results can find applications in foundations of quantum theory, complexity studies 
and classical simulations of quantum systems. 

PACS numbers: 03.65.Ta, 03.65.Ud, 03.67.Lx, 89.70.Eg 



In classical physics, the position and momentum of a 
particle determine the outcomes of all possible measure- 
ments that can be performed upon it. They define a 
deterministic classical state. If the state is not fully ac- 
cessible, a general probabilistic classical state is a mixture 
of the deterministic states, arising from the inaccessibil- 
ity. Since quantum mechanics gives only probabilistic 
predictions, it was puzzling already to the fathers of the 
theory whether it can be completed with an underlying 
classical-like model jT] . The quantum probabilities would 
then arise from an inaccessibility of some hidden variables 
(HV) describing analogs of deterministic classical states, 
the hidden-variable states, which determine the results 
of all quantum measurements. 

Since the seminal work of Kochen and Specker (KS), it 
has been known that HV models must be contextual [SJ. 
On the operational level, the contextual HV models can- 
not be distinguished from quantum mechanics. However, 
one may ask how plausible these models are in terms of 
resources, e.g., how many HV states (also called the "on- 
tic states" [3]|I1[5]) they require. In addition to the fun- 
damental question of the minimal HV model for a quan- 
tum system, this research is motivated by problems in 
quantum information theory. In particular, HV models 
allow a fair comparison between complexities of quantum 
and classical algorithms [5J [7] , as a quantum algorithm 
can now be represented by a classical circuit. 

For an infinite number of measurement settings, al- 
ready a single qubit requires infinitely many HV states, 
the result proved by Hardy [Bj and, in a different con- 
text, by Montina [9], [10]. However, these authors did not 
consider the scaling of the number of HV states with the 
number of measurements. Harrigan and Rudolph found 
a deterministic HV model that requires exponentially 
many HV states to simulate results of the finite set of 
measurements on all quantum states |llj . Our construc- 
tion also provides such models and consumes at most a 
polynomial number of HV states, bringing exponential 
improvement. In the limit of infinitely many measure- 



ments, the number of HV states for an indctcrministic 
model scales linearly with the number of measurements. 
Moreover, the number of real parameters that specify 
these HV states saturates the lower bound derived by 
Montina |10j and, consequently, is the minimal number 
possible. Our method also allows a universal generaliza- 
tion of the Spekkens model [5] . 

Consider a finite number, N, of projective measure- 
ments on a d-level quantum system in a state p. The 
probability to observe the rth result in the nth measure- 
ment is Pr(p) — Tr[pflr ], where fli n ^ is a projector 
on the rth orthogonal state of the nth measurement, i.e., 
r = 1, ...,d and n = 1, N. We form a (/-dimensional 
vector, pM = (jp^ , . . . ,p)^) T , composed of the prob- 
abilities for distinct outcomes in the nth measurement. 
For the set of measurements, we build a dTV-dimensional 
preparation vector, p = (p^, ...,p^ N,> ) T [T2]. The deter- 
ministic HV states predetermine the results of all mea- 
surements and can be represented as a diV-dimensional 
vector 

O ri ... rN = (0,...,1,...,0|...|0,...,1,...,0) T , (1) 

where r n is the position of 1 in the nth sequence (r n = 
0, . . . , d — 1 indicates that outcome r n occurs in the nth 
measurement). The space of all HV states, A, is formed 
by classical mixtures of d N deterministic states O ri ...r N - 
A set of k quantum states p\,...,p K has a HV 
model for N measurements, if one can find L vectors 
Ox,..., Or, g A such that 

L 

p(Pfe) = ^aii^Oi, for all k = 1, k (2) 

2 = 1 

where cv/(fc) > and J2i a i(k) = 1- The model is 
called deterministic if all O/ are deterministic HV states; 
otherwise, it is called indeterministic. The model is 
preparation-universal, if the HV states simulate any phys- 
ical state p, and it is measurement-universal if they sim- 
ulate any measurement. 
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Formally, the set A is a convex polytope in Hl dN having 
the states O; as vertices. Since all probabilities satisfy 
» 



show a nonuniversal model, only for the eigenstates 
of these measurements: ±m x , ±m y , ±m z . The cor- 

(+) 



< Pr < 1, any preparation vector p(p) lies inside this responding preparation vectors are: pi ' = (1, |, 1) 
polytope and has a HV model. We study the number of 
HV states required for the model. 

We begin with a specific deterministic HV model for 
a two-level quantum system (qubit) which we shall often 
refer to later on. An arbitrary state of a qubit can be rep- 
resented as p = 5(1+^1=1 x i&i)i where fri's are the Pauli 
matrices and x = (x\, X2, X3) T is a Bloch vector, in a unit 
ball |x| < 1. A set of N projective measurements, with 
2N outcomes (states on which the qubit is projected), 
is described by 2N unit vectors ±ixii, ■ • ■ , ±mjv on the 
Bloch sphere. The preparation vector for these directions 

' l±mix 



pi > = (0,i,i), p<, +) = (1,1, i), pi _) = (|,0, i), and 

pi +) = (§, 1, 1), pi _) = (|, |,0). Applying the method 
of ([3]) to each of these preparation vectors, one finds 
that L = 4 HV states are sufficient for the simulation: 
O - (1,1, 1) T , Ch = (1,0,0) T , 2 = (0,1, 0) T , and 
O3 = (0,0, 1) T . These four states, together with their 
decomposition of the preparation vectors, 



is p(x) = 
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). Since the probability 
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for the measurement — m is fully determined by the one 
for the +m, one can reduce ("compress") preparation 
vector to p(x) = (i±=U5, . . . , ±±™4ES). Similarly, the 
deterministic HV states are reduced to A^-dimensional 
vectors O ri ... rjv = (n, . . . , r]y) T , where r„ = 0,1. The 
(reduced) space A is a hypercube in N dimensions, with 
2 N vertices defined by these states. By Caratheodory's 
theorem [15] for each vector p(x) = (pi, . . . ,pn) T , one 
can identify N + 1 HV states the convex hull of which 
contains p(x). For a given x, the vector p(x) can be writ- 
ten as a permutation of a reordered preparation vector 
p^(x) wherein the probabilities appear in increasing or- 
der, p[ < p\ < • • ■ < and the latter can be expressed 
in terms of N + 1 HV states as 
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where the columns of the displayed matrix are the HV 
states. The expansion coefficients are 



a 



= 1 - 

= P l n 



-pi-i 



for 



■A, 

n = 2,...,N, 



(4) 



and, due to the ordering of probabilities, the coefficients 
are all positive and sum up to 1. One can suitably per- 
mute the rows in matrix given by ^ to bring the proba- 
bilities in order given by p(x). Thus, p(x) can be written 
as a convex combination of A^+l columns (HV states) of a 
reordered matrix. The number of A^+l states can be fur- 
ther reduced. E.g., for two equal probabilities, pi = p 2 , 
the number of HV states is decreased because a-i = 0. If, 
say, P2 — 1 — Pi, one can exchange mi — > —mi, such that 
the probabilities become equal, leading to another reduc- 
tion. Importantly, different quantum states are generally 
modeled by different sets of N + 1 HV states. 

As an illustrative example, consider a model for three 
complementary measurements along m x ,mj,,m z . We 
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2° 1 + 2° 3 < 
2° 1 + 2° 2 - 



(5) 



are equivalent to the toy model of Spekkens [5]. 

We give a constructive proof that a preparation- 
universal simulation of N quantum measurements on a 
qubit can be achieved with the number of HV states 
that is polynomial in N. Let Ai denote a polytope 
formed as a convex hull of the measurement settings, 
Ai = conv{±mi, . . . , ±mjv}- Its dual polytope is a set 
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l<m„y<l,n=l...iV}. (6) 



The polytope Ai lies inside the Bloch sphere and its 
dual contains the sphere. Therefore, every Bloch vec- 
tor can be written as a convex combination of the ver- 
tices, y;, of the dual polytope, x = X^ a z( x )yz- The 
components of the measurement vector can now be de- 
composed as p„(x) = J2i a i( x )i(l + m nYz)- Accord- 
ing to the definition of the dual polytope, the quantity 
1(1 + m n y;) S [0,1] and can be interpreted as the nth 
component (probability) of the Zth HV state. Since the 
Bloch vectors corresponding to projections onto orthog- 
onal states sum up to the zero vector, the correspond- 
ing probabilities assigned by a HV state sum up to 1, 
as it should be. Thus, the set of HV states correspond- 
ing to vertices of the dual polytope is sufficient for a 
preparation-universal HV model. Note that this model 
can in general be indeterministic. In such a case, each 
indeterministic HV state can be further reduced into at 
most N — 2 deterministic HV states, according to pj. 
The reason for N — 2, and not N + 1, states stems 
from the observation that a vertex of the dual polytope 
saturates at least three of the inequalities defining the 
polytope (at least three facets have to meet at each ver- 
tex), i.e., the corresponding probability is 1 or 0, and 
reduces the number of required deterministic HV states. 
Finally, the total number of HV states required for an 
indeterministic model is L < F, and for a deterministic 
model is L < (N — 2)F, where F is the number of ver- 
tices of the dual polytope or, equivalently, the number of 
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facets of the measurement polytope. A convex polytope 
with 2N vertices (in three-dimensional space) can have 
N + 2 < F < 4(iV - 1) facets [13 , which implies that 
indeterministic HV models require at most a number of 
HV states that is linear in N, and deterministic ones 
require quadratic number of HV states. 

Using the dual polytope approach, we generalize 
Spekkens' model [5], originally formulated to explain 
the measurement results on the eigenstates of the three 
complementary directions, to the preparation-universal 
model. For these directions, the measurement polytope 
is an octahedron, see Fig. [TJa). The dual polytope is a 
cube, whose interior forms the whole space of HV states, 
with the vertices being the deterministic states. Another 
interesting example is illustrated in Fig.JlJb). 

The dual polytope approach can be applied to arbi- 
trary preparation vectors. However, efficient simulations 
are only expected for highly symmetric polytopes. For 
this reason, we move to more complicated Platonic solids 
and general symmetry considerations. 

Consider a set of measurement directions ±mi, . . . , 
im^r, which is generated by a group; e.g., an octahe- 
dron and a cube can be generated via the chiral octa- 
hedral group O with 24 rotations. Generally, if G is a 
symmetry of the measurement polytope, A4, it is also a 
symmetry of its dual, T>w, i.e., the dual polytope can 
also be generated by G. The group action permutes the 
vectors ±m„ as well as vertices of the dual polytope. 
Since the last are related to the HV states, we can define 
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FIG. 1: Preparation-universal HV models and dual poly- 
topes. (a) The vertices of the octahedron inside the Bloch 
sphere define the three complementary qubit measurements. 
A preparation-universal HV model for these measurements 
requires eight HV states, which are written near their repre- 
sentative vertices of the cube containing the sphere. It gen- 
eralizes the Spekkens model [5], which is not universal and 
utilizes only four out of eight states (see main text). Their 
corresponding vertices span a tetrahedron inside the cube, 
which does not contain the whole Bloch sphere, (b) Here, 
the measurement directions form a cube inside the sphere. Al- 
though more measurements are to be simulated, the universal 
HV model requires only six HV states, which are written near 
their representative vertices of the octahedron containing the 
sphere. 



the permutation representation of the group in the HV 
space, Dp(G). The HV state, h(y'), corresponding to a 
vertex of a dual polytope, y' = gy, which is generated by 
g G G acting on an initial vertex, y, can be found using 
the group representation: 

h(.gy) = Dp(s)h(y). ( 7 ) 

Decomposing h(y) into deterministic HV states brings 
to the form h(<?y) = £ J=1 a,£> P ( ff )0,. There- 
fore, the set of deterministic HV states required for the 
preparation-universal model is the union of a number of 
group orbits {Dp(g)Oi\g e G}. Because of the symme- 
tries involved, the minimal number of HV states cannot 
be smaller than the number of elements in the smallest 
orbit. 

Let us consider two other Platonic solids, the icosahe- 
dron and the dodecahedron [17]. Both of them posses the 
same symmetry, the chiral icosahedral group X, with 60 
rotations. Consider the icosahedron as the measurement 
polytope, N = 6. Its dual, the dodecahedron, has 20 ver- 
tices corresponding to indeterministic HV states that can 
be further reduced to deterministic HV states. The total 
number of possible deterministic HV states is 2 6 = 64 in 
this case. We have found four different orbits of action of 
X with 12, 12, 20, 20 different elements, respectively. Only 
one orbit, with 20 elements, gives deterministic states for 
universal simulation. For N — 10 measurement settings, 
the dodecahedron is the measurement polytope. Its dual, 
the icosahedron, has 12 vertices. The total number of 
possible deterministic HV states is 2 10 = 1024, which is 
partitioned into 24 different orbits: 2 with 12 elements, 
8 with 20, and 14 with 60 elements. The two lowest 
orbits are suitable for the universal model. Thus, the 
minimal deterministic model, among all HV models ob- 
tained through the dual polytope construction, requires 
only 24 HV states, twice the number of vertices of the 
dual polytope. 

The presentation so far was limited to qubits. How- 
ever, a similar line of reasoning applies to any d- level 
quantum system. In the general case, Pauli operators 
have to be replaced by generalized Gell-Mann operators, 
Xi, which naturally leads to the generalized, D = d 2 — 1 
dimensional, Bloch representation. An arbitrary quan- 
tum state, p — t[1 + {d — 1) Yl^L^ is now repre- 
sented by a generalized Bloch vector, x, with components 
Xi = Tr(/5Ai). We normalize the Gell-Mann operators 
as Tr(AiAj) = jzj<5ij, such that pure quantum states 
are represented by normalized generalized Bloch vectors. 
Contrary to the qubit case, not every unit vector corre- 
sponds to a physical state. The probability of an out- 
come associated with a projector on a state represented 
by m„, in a measurement on a state represented by x, is 
p n (x) = h[l + (d— l)(m n x)]. The requirement of positive 
probabilities reveals that, e.g., the vector x = m n does 
not represent a physical state. 
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In analogy to the dual poly tope, for a set of dN 
preparaion vectors, representing TV d-valued observables, 
we introduce a convex polytope the interior of which in- 
cludes all vectors y leading to physically allowed proba- 
bilities Pn (y) £ [0,1]: 

V M = {y e R D \ - < m„y < 1, n = 1, ...,dN}. (8) 

Among others, this polytope contains all the vectors of 
quantum states. The generalized Bloch vectors corre- 
sponding to a complete set of orthogonal quantum states 
sum up to the zero vector, implying the probabilities as- 
signed by a HV state for different outcomes of any mea- 
surement sum up to 1, as it should be. Again, the vectors 
of quantum states can be expressed as a convex combi- 
nation of vertices of Vmi an d their number gives the 
upper bound on the amount of HV states sufficient for 
preparation-universal simulation. The polytope Vm 1& 
specified by q = 2dN linear inequalities, two inequalities 
for each vector m n , and its maximal number of vertices 
is given by L < (*Z S D ) + where 6 = [(D + l)/2j, 

5' = [(D + 2)/2j, and [a; J is the integer part of x [13]. In 
the special case of a qubit, the dual polytope is defined 
by 2N, and not 4JV, inequalities because the two bounds 
of Eq. ^ are the same for the vectors ±m n . Since the 
binomial coefficient increases with a, L < 2( 9 ~^,). 

Using (I) = ( a ° fc ), we have L < 2(^), and since 
(?) < a b /bl, the maximal number of vertices is polyno- 
mial in N, L ~ (2dN - 5) D - 5 . The related HV states 
can in general be indeterministic, and each of them can 
be decomposed to O(N) deterministic HV states, using 
decomposition ^ in the dN dimensional space A. There- 
fore, for any system, the number of (in) deterministic HV 
states required for a preparation-universal simulation is 
polynomial in N. 

In the limit of infinitely many measurements, our 
method gives (preparation and measurement) universal 
models with the minimal number of HV states. As proved 
by Montina, in this limit the optimal model requires 
2{d — 1) real parameters to describe the HV states [10] , 
We show that for an infinite number of settings the set of 
universal HV states converges to the set of pure quantum 
states, which is known to be parameterized by 2(d — 1) 
real numbers. First, consider a finite set of projectors H n 
with n = l,...,dN, and the corresponding polytope ^ 
in the Hilbert-Schmidt space of Hermitian operators with 
unit trace. The operators of its vertices, yi, correspond 
to the HV states, i.e., for all n, Tr(yjII n ) gives the prob- 
ability that is assigned by the HV state, of the outcome 
associated with projector fl„. For other projectors, not 
within the set of dN, the trace does not have to represent 
a probability and therefore the set of operators yi is larger 
than the set of quantum states [18]. However, in the limit 
of infinitely many measurements, Tr(y;II n ) £ [0, 1] for all 
possible projectors; therefore, the eigenvalues of yi's lie 
within the [0, 1] interval. Since Tr(y;) = 1, the operators 



pi are just quantum states and the HV states correspond- 
ing to pure quantum states are universal. Their number 
scales linearly with N, because N measurements corre- 
spond to dN projectors and each of them represents one 
HV state (and also one pure quantum state). 

Regarding the polytope Vm m the space of Hermi- 
tian operators allows for an easy generalization of our ap- 
proach to POVM measurements. POVM elements, E n , 
are positive operators being vertices of a measurement 
polytope. The polytope Vm includes all the unit-trace 
operators y for which Tr (yE n ) £ [0,1]. Since for all quan- 
tum states Tr(pE n ) £ [0,1], the polytope Vm contains 
all of them and, as before, its vertices define HV states. 

For a d-level system the KS argument disqualifies non- 
contextual HV theories [2], and one might wonder how 
contextuality enters our models. Consider the KS ar- 
gument of Peres [14] . It involves 33 different vectors in 
M 3 , which belong to 16 different orthogonal triads. Non- 
contextuality requires a value associated with a single 
vector to be the same irrespectively of other vectors in the 
triad. In the present models, the results of 16 different 
measurements are described by HV states with 3-16 = 48 
components; i.e., a value assigned to the same vector can 
depend on the other vectors in the triad. 

In conclusion, we proved that a preparaion-universal 
HV model of the results of N quantum measurements re- 
quires at most a number of HV states which is polynomial 
in A. In the limit of infinitely many measurements, our 
method gives optimal preparation- and measurement- 
universal HV models, with the minimal number of real 
parameters describing the HV states. There is no HV 
model that would require less HV states than the model 
in which every quantum state is associated with a HV 
state [10] . Furthermore, since there are infinitely many 
measurements that can be performed on a quantum sys- 
tem, its HV description requires infinitely many HV 
states. This "ontological baggage" [5] can be seen as 
an argument against the HV approach because it is ex- 
tremely resource demanding already for a single qubit. 
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